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2 ( 57 ) Abstract: Improved Wavefront Coding Optics (106, 601, 602), which apply a phase profile to the wavefront of light from an 
object to be imaged, retain their insensitivity to focus related aberration, while increasing the heights of the resulting MTFs and 

^ reducing the noise in the final images. Such improved Wavefront Coding Optics have the characteristic that the central portion of 
the applied phase profile is essentially flat (or constant), while a peripheral region of the phase profile around the central region 

^ alternately has positive and negative phase regions relative to the central region. 
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MTF- IMROVED OPTICAL SYSTEM EMPLOYING PHASE MASK WITH UNCHANGED PHASE IN 
CENTRAL REGION 

PACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION: 

This invention relates to improved Wavefront Coding Optics for controlling 
focus related aberrations, and methods for designing such Wavefront 
Coding Optics. 

DESCRIPTION OF THE PRIOR ART: 

Wavefront Coding is a relatively new technique that is used to reduce the 
effects of misfocus in sampled imaging systems through the use of 
Wavefront Coding Optics which operate by applying aspheric phase 
variations to wavefronts of light from the object being imaged. Image 
processing of the resulting images is required in order to remove the spatial 
effects of the Wavefront Coding. The processed images are sharp and 
clear, as well as being relatively insensitive to the distance between the 
object and the detector. Wavefront Coding is also used to control general 
focus related aberrations to enable simplified design of imaging systems as 
well as to provide anti-aliasing in sampled imaging systems. 

The Wavefront Coding Optics taught and described in prior art, beginning 
with U.S. Patent No. 5,748,371, issued May 5, 1998, were discovered by trial 
and error. The first operative Wavefront Coding mask applied a cubic 
phase function to the wavefront from the object. It was known that 
Wavefront Coding Optics, like the cubic mask, needed to apply aspheric, 
asymmetric phase variations to the wavefronts. 

Prior art related to Wavefront Coding systems includes a fundamental 
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description of Wavefront Coding (U.S. Patent No. 5,748,371), description of 
Wavefront Coding used for anti-aliasing (Anti-aliasing apparatus and 
methods for optical imaging, U.S. Patent #6,021,005, February 1, 2000), 
use of Wavefront Coding in projection systems (Apparatus and methods for 
5 extending depth of field in image projection systems, U.S. Patent 

#6,069,738, May 30, 2000), and the combination of Wavefront Coding and 
amplitude apodizers (Apparatus and method for reducing imaging errors in 
imaging systems having an extended depth of field, U.S. Patent #6,097,856, 
August 1, 2000). 

1 0 The layout of a conventional Wavefront Coded imaging system is shown in 
Figure 1. Imaging Optics 104 collects light reflected or transmitted from 
Object 102. Wavefront Coding Optics 106 modify the phase of the light 
before detector 108. Wavefront Coding Optics 106 comprise a cubic mask. 
Detector 108 can be analog film, CCD or CMOS detectors, etc. The image 

1 5 from detector 108 is spatially blurred because of Wavefront Coding Optics 

106. Image processing 110 is used to remove the spatial blur resulting in a 
final image, i.e., image processing 110 removes the Wavefront Coding 
applied by optics 106, thereby reversing the effects of optics 106, other 
than the increase in depth of field and depth of focus. The image before 

2 0 and after Image Processing 110 also is very insensitive to misfocus 

aberrations. These misfocus aberrations can be due to the Object 102 being 
beyond the depth of field of the Imaging Optics 104, the detector 108 being 
beyond the depth of focus of the Imaging Optics 104, or from Imaging 
Optics 104 having some combination of misfocus aberrations such as 
25 spherical aberration, chromatic aberration, petzval curvature, astigmatism, 
temperature or pressure related misfocus. 

Figure 2 describes a rectangularly separable prior art Wavefront Coding 
phase tunction that produces an extended depth of field. This phase 
function is a simple cubic phase function that is mathematically described, in 
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normalized coordinates, as: 

cubic-phase(x.y) = 12 [ x a + y 3 ] 
Ixl < 1, lyl < 1 

Other related forms of the cubic mask are described as: 

5 cubic-related-forms(x,y) = a[sign{x)lxl b + sign(y)lyl b ], 

where 

Ixl £ 1, lyl £ 1, 

and 

sign(x) = +1 for x £ 0, sign(x) = -1 otherwise 

1 0 These related forms trace out "cubic like" profiles of increasing slopes near 
the end of the aperture. 

The top plot of Figure 2 describes a 1D slice along an orthogonal axis of the 
cubic phase function. The lower plot of Figure 2 describes the contours of 
constant phase of this cubic phase function. 

1 5 Figure 3 shows MTFs as a function of misfocus for a system with no 

Wavefront Coding and for a system with the conventional Wavefront Coding 
cubic phase function of Figure 2. The normalized misfocus values are the 
same for both systems and are given as v = { 0 , 2 , 4}, where \j/ = [ 2 pi 
W 20 ], and where W 20 is the conventional misfocus aberration coefficient in 

20 waves. MTFs with no Wavefront Coding (302) are seen to have a large 
change with misfocus. MTFs with the rectangularly separable cubic phase 
function (304) are seen to change much less with misfocus then the system 
with no Wavefront Coding. 

A non-separable prior art form of Wavefront Coding Optics, in normalized 
25 coordinates, is: 

non-separable-cubic-phase( p, q ) = p 3 cos( 3 q ) 
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Ipl < 1, 0 < q < 2 pi 

This phase function has been shown to be useful for controlling misfocus and 
for minimizing optical power in high spatial frequencies, or antialiasing. When 
using a digital detector such as a CCD or CMOS device to capture image 
5 108, optical power that is beyond the spatial frequency limit of the detector 

masquerades or "aliases" as low spatial frequency power. For example, say 
that the normalized spatial frequency limit of a digital detector is 0.5. As 
seen from Figure 3, the in-focus MTF from the conventional system with no 
Wavefront Coding can produce a considerable amount of optical power 
1 0 beyond this spatial frequency limit that can be aliased. By adding misfocus 
to the system without Wavefront Coding the amount of high spatial 
frequency optical power can be decreased, and aliasing reduced, as is well 
known. When using conventional Wavefront Coding, as shown in Figure 3, the 
amount of optical power that can be aliased can be decreased (304) 

1 5 compared to the system without Wavefront Coding (302). 

Image Processing function 110 essentially applies amplification and phase 
correction as a function of spatial frequency to restore the MTFs before 
processing to the in-focus MTF from the conventional system with no 
Wavefront Coding after processing, or to some other application specific 
20 MTF, if desired. In effect, the Image Processing function of Figure 1 removes 

the Wavefront Coding blur in the detected image. 

In practice the amplification applied by the Image Processing function 
increases the power of the deterministic image but also increases the power 
of the additive random noise as well. If Image Processing 110 is 

2 5 implemented as a linear digital filter then a useful measure of the increase of 

power of the additive random noise is called the Noise Gain of the digital 
filter. The concept of "noise gain" is commonly used in radar systems to 
describe the amount of noise power at the output of radar digital 
processors. Nonlinear implementations of image Processing 110 have 
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similar types of noise-related measures. The Noise Gain for a digital filter is 
defined as the ratio of the root-mean-square (RMS) value of the noise after 
filtering to the RMS value of the noise before filtering. In general the Noise 
Gain is nearly always greater than one in Wavefront Coded systems. 
5 Assuming that the additive noise is uncorrected white gaussian noise, the 

Noise Gain of a two dimensional linear digital filter can be shown to be equal 
to: 

Noise Gain = sqrt[ Z Z f(i,k) 2 ] = sqrt[ Z Z IFtw^w^i 2 ] 

where 

10 I ZZ f(i,k) ] = F(0,0) = 1.0, 

f(i,k) is a spatial domain digital filter, F(w f ,w k ) is the equivalent frequency 
domain digital filter, and the first sum is over the index i or k and the 
second sum is over the other index. Indices (i,k) denote spatial domain 
coordinates while indices (Wj,w k ) denote frequency domain coordinates. The 

1 5 constraints that the sum of all values of the filter and the zero spatial 

frequency filter value both equal unity ensures that the zero spatial 
frequency components of the image (the background for example) are 
unchanged by the image processing. 

Wavefront Coded MTFs that have the highest values require the least 
20 amplification by the digital filter and hence the smallest Noise Gain. In 

practice the Wavefront Coding Optics that produce MTFs that have small 
changes over a desired amount of misfocus and also have the highest MTFs 
are considered the best and the most practical Optics for Wavefront 
Coding. Optics that produce MTFs that have small changes with misfocus 

2 5 but also very low MTFs are impractical due to very large Noise Gain of the 

resulting digital filters. Digital filters with large Noise Gain will produce final 
images that have unnecessarily high levels of noise. 
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While the conventional cubic Wavefront Coding mask does operate to 
increase depth of field and control focus related aberrations, there remains 
a need in the art for improved Wavefront Codings Optics, which retain the 
capacity to reduce focus-related aberrations, while also producing high value 
5 MTFs. There also remains a need in the art for methods of designing such 

improved Wavefront Coding Optics. 

SUMMARY OF THE INVENTION 

An object of the present invention is to provide improved Wavefront Coding 
Optics, which retain the capacity to reduce focus-related aberrations, while 
1 0 producing high value MTFs, and to provide methods of designing such 

improved Wavefront Coding Optics. 

Improved Wavefront Coding Optics according to the present invention, which 
apply a phase profile to the wavefront of light from an object to be imaged, 
retain their insensitivity to focus related aberrations, while increasing the 
1 5 resulting MTFs. Such improved Wavefront Coding Optics have the 

characteristic that the central portion of the applied phase profile is 
essentially constant, while the edges of the phase profile have alternating 
negative and positive phase regions (e.g. turn up and down at respective 
profile ends). 

20 To achieve higher MTFs, control misfocus and misfocus aberrations, and 

improve antialiasing characteristics, the central group of rays should be left 
unmodified by the Wavefront Coding Optics. In order to increase the light 
gathering (and possibly spatial resolution) of the full aperture system, the 
outer rays need to be modified. Only these outer rays need to be modified 

25 in order to increase the light gathering while keeping the depth of field 

and/or aliasing characteristics constant. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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Figure 1 (prior art) shows a conventional prior art Wavefront Coding 
imaging system. 

Figure 2 (prior art) shows a 1D plot of the Wavefront Coding cubic phase 
function of Figure 1 and a contour plot of the 2D representation of this 
5 function. 

Figure 3 (prior art) shows MTFs as a function of misfocus for a system with 
no Wavefront Coding and with the conventional Wavefront Coding cubic 
phase function of Figure 2. 

Figure 4 (prior art) shows the ray paths for a system with no Wavefront 
1 0 Coding and with the conventional Wavefront Coding cubic phase function 

shown in Figure 2. 

Figure 5A shows the ray paths for a linear region phase function according 
to the present invention. 

Figure 5B shows the ray paths for a cubic region phase function according 

1 5 to the present invention. 

Figure 6 shows 1D plots of the linear region and cubic region phase profiles 
of Figures 5a and 5b, as well as contour plots of the 2D representations of 
these Wavefront Coding functions. 

Figure 7 shows MTFs as a function of misfocus for a system with no 

2 0 Wavefront Coding and with the linear region and cubic region Optics of 

Figures 5a and 5b. 

Figure 8 shows a contour plot of an improved non-separable phase function 
according to the present invention, a plot of a 1D slice through the non- 
separable phase function, and the MTFs as a function of misfocus for a non- 
2 5 separable cubic sector phase function. 



7 



WO 03/021333 PCT/US02/09702 

Figure 9 (prior art) shows a drawing of a conventional Cooke triplet lens. 

Figures 10A and 10B (prior art) show contour plots of the exit pupils and 
corresponding MTFs as a function of field angle for the full aperture (1 OA) 
and stopped down (10B) Cooke triplet lens of Figure 9. 

5 Figures 11A and 11B show contour plots of the exit pupils and 

corresponding MTFs as a function of field angle for the prior art cubic phase 
Wavefront Coding system of Figure 1 (11 A) and the improved cubic region 
Wavefront Coding system of Figure 8 (11B). 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

1 0 There are an infinite number of Wavefront Coding Optics that will reduce the 

variation in the resulting MTFs and PSFs of a given optical system due to 
misfocus or misfocus aberrations. Many of the possible optics are 
impractical in that the required Image Processing function 110 used to 
remove the Wavefront Coding blur from detected images would amplify the 

1 5 additive noise in practical images beyond an acceptable level. Improved 

forms of Wavefront Coding Optics that can control misfocus and misfocus 
aberrations, that can lead to higher MTFs, and that have improved 
antialiasing characteristics, as well as new methods of Wavefront Coding 
design, are shown in Figures 5 through 9. Use of these improved forms of 

20 Wavefront Coding Optics and design methods to control the misfocus 
aberrations with a Cooke triplet lens are shown in Figures 9 through 11. 

The improved Wavefront Coding Optics according to the present invention 
share the characteristic that the central region of the applied phase profile 
is essentially constant, while the edges of the phase profile have alternating 
25 positive and negative phase regions. Such Wavefront Coding Optics retain 
the ability to reduce focus related aberrations, but also have considerably 
higher MTFs than traditional Wavefront Coding Optics, thus reducing noise in 
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the final images produced. 

Wavefront Coding Optics (in the form of aspheric optical elements) are 
placed at or near the aperture stop of optical systems (or at or near an 
image of the aperture stop) in order to redirect light rays as a function of 
5 spatial position. The aspheric optical elements can be constructed of 

optical glass or plastic with varying thickness and/or index of refraction. 
The optics can also be implemented with shaped mirrors, spatial light 
modulators, holograms, or micro mirror devices. U.S. Patent No. 
6,021,005, issued February 1, 2000 and entitled "Antialiasing Apparatus and 
1 0 Methods for Optical Imaging," provides descriptions of a variety of devices 

for applying variations to a wavefront of light from an object. 

Light rays from an ideal thin lens with no Wavefront Coding converging to 
focus 50mm from the lens are shown on the top graph of Figure 4. All 
rays from the ideal system without Wavefront Coding travel towards the 

1 5 best focus position on the optical axis. The light rays from a conventional 

(prior art) rectangularly separable cubic phase system are shown on the 
bottom graph of Figure 4. Notice that the rays from the top half of this 
lens cross the optical axis beyond the best focus point of the conventional 
lens (or 50mm). The rays from the bottom half of this lens cross the 

2 0 optical axis before the best focus point of the conventional lens. 

Since no two rays from the prior art Wavefront Coding system cross the 
optical axis at the same point, every ray of the prior art Wavefront Coding 
cubic phase system is being modified when compared to the system not 
using Wavefront Coding (except the on-axis zero-slope ray). To achieve 
25 higher MTFs, control misfocus and misfocus aberrations, and improve 

antialiasing characteristics, the central group of rays should be left 
unmodified. 

Consider a full aperture and a stopped down conventional optical system 
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that does not use Wavefront Coding. Assume that a particular application 
has a depth of field (or depth of focus) and/or antialiasing requirement that 
the full aperture system cannot meet. It is well known that stopping down 
the lens will increase the depth of field of the system while reducing the 
5 potential spatial resolution of the optics. Stopping down the aperture also 

reduces the optical power available at the detector. There is a particular 
stopped down aperture where the imaging system has a best fit to the 
depth of field and/or antialiasing characteristics required. The light rays 
that pass through the aperture of the stopped down system are then 
1 0 considered suitable from a depth of field and/or antialiasing perspective for 
the particular application. 

In order to increase the light gathering (and possibly spatial resolution) of 
the full aperture system, the rays of the full aperture system that lay 
outside the stopped down aperture need to be modified. Only those rays 

1 5 that lay outside the stopped down aperture need to be modified in order to 
increase the light gathering while keeping the depth of field and/or aliasing 
characteristics constant. When the rays within the stopped down aperture 
are modified, as in all prior art Wavefront Coding Optics, the resulting MTFs 
are not as high as possible, the resulting Noise Gains are not as low as 

20 possible, and the resulting images are more noisy than necessary. 

By not modifying the central rays of Wavefront Coding systems the resulting 
MTFs below the spatial frequency limit of the digital detector can be 
increased over prior art Wavefront Coding systems. Only the rays outside 
of the central region of the aperture need to be modified in order to control 
25 misfocus or misfocus aberrations. The central region is defined as the 

general region of a stopped down aperture where a system would have a 
suitable depth of field, depth of focus, or antialiasing characteristics for a 
particular application, albeit with reduced light gathering and spatial 
resolution. 
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Improved rectangularly separable Wavefront Coding Optics can be most 
generally described mathematically through phase functions, in normalized 
coordinates, that have the form: 

phase(x,y) = I[ U(lxl/A xi ) G xi (x) + U(lyl/A yl ) G yj (y) ] 

vyhere 

Ixl £ 1, lyl s 1 
1 = 1,2 N 

and where 

U(z) = 1 if z >1, U(z) = 0 otherwise 
0< A^l, 0< A yi <1 

The sum is over the index i. The function U(lxl/A x ) is a zero/one step 
function that has a value of zero inside of an aperture of length 2 A x and 
has a value of one outside of this aperture. The functions G x and G y are 
general functions that modify the rays of the system outside of the specific 
aperture defined by A x and A y . The aperture shape in this form is described 
by a rectangle for mathematical convenience, but in general can be 
described by any closed shape. For example, instead of a rectangular 
aperture, a circular, elliptical, or multi-sided polygonal aperture can also be 
be used. 

In order for the phase function to control misfocus effects, the phase 

functions G v and G„ should be designed so that groups of rays from specific 
* y 

regions of the aperture cross the optical axis either before or after the 
point of best focus when no Wavefront Coding is used. 

With these concepts, a number of improved rectangularly separable 
Wavefront Coding Optics can be formed. These optics differ in the 
composition of the general functions G x and G y . For example, a linear phase 
region optical system can be described as: 

1 1 
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linear-phase-region(x,y) = U(lxl/A x ) sign(x)(lxl«A x )/(1-A x ) 
+ U(lyl/A y ) sign(y)(lyl-A y )/(1-A y ) 

where 

Ixl <> 1, lyl < 1 

5 U(z) = 1 if z <1, U(z) = 0 otherwise 

sign(z) = +1 for z £ 0, sign(z) = -1 otherwise 
0< A x <1, 0< A y <1 

The linear phase region system has zero phase inside of the aperture 
defined by A x and A y with linearly changing phase as a function of spatial 
1 0 position variables x and y outside of this aperture. The linear phase region 
system delivers extended depth of field with high MTFs below the spatial 
frequency limit of the digital detector. This system also has a very simple 
physical form that can be preferred over smoothly varying forms because 
of the fabrication process or physical implementation being used. 

1 5 For example, if fabrication of the Wavefront Coding Optics is done with a 

precision milling machine, the linear phase region optics would be preferred 
over the cubic phase optics because the linear phase region optics have only 
fixed surface slopes while the cubic phase optics have continuously changing 
surface slopes. Often, optical elements with fixed surface slopes are easier 
20 to fabricate than elements with continuously changing slopes. If the 

Wavefront Coding optical surface is implemented with micro mirrors or 
similar, the small number of surface slopes of the linear phase region system 
can be more easily implemented than the continuously changing surface 
slopes of the prior art cubic phase optics. 

2 5 Another version of improved rectangularly separable Wavefront Coding 

Optics is mathematically described by the sum of powers region phase 
function. This phase function is defined as: 
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sum-of-powers-region(x,y) = 

Z U(lxl/A xi ) ttj sign(x) [ (Ixl -A xi )/(1-A xi ) ] Sj 

+ Z U(lyl/A yl ) Xi sign(y) [ (lyl -A yj )/(1-A yl ) ] 6i 

where 

5 Ixl < i, lyl < 1 

i = 1,2, ... . N 
U(z) = 1 if z < 1 , U(z) = 0 otherwise 
sign(z) = +1 for z £ 0, sign(z) = -1 otherwise 
0< A xl <1, 0< Ay { <1 

-j o and where the summations are over the index i. 

Examples of the linear phase region systems and the sum of powers region 
systems are found in Figures 5 t 6, and 7. The graph of Figure 5A shows the 
ray paths from one dimension of the linear phase region system. The graph 
of Figure 5B shows the ray paths for one dimension of a cubic phase region 

1 5 system. A cubic phase region system is sum of powers phase region 

system with only a single term and the exponent on this term having the 
value of 3, or Bj = b } = 3, and <Xj = Xj = 0 for i * 1 . 

From Figure 5A we see that the rays from the central region of the lens 
(with this central region extending from -2.5 to +2.5mm) are not modified 
20 by the linear phase region system. All the unmodified rays travel towards 
the optical axis at the best focus position 50mm to the right of the lens. 
Only the outer rays from the lens, those beyond +/-2.5mm, are modified by 
the linear phase region function. The rays from the top of the lens are 
modified so that they cross the optical axis at points beyond the best focus 

2 5 position. The rays from the bottom of the lens cross the optical axis at 

points before the best focus point. 

From Figure 5B, the cubic phase region system also does not modify the 
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central region (extending from -2.5mm to +2.5mm) of the lens. The 
remaining rays are modified in a spatially varying manner such that one set 
of rays (from the top of the lens) crosses the optical axis at points beyond 
the best focus point and another set of rays (from the bottom of the lens) 
5 crosses the optical axis at points before the best focus position. 

The linear phase region system and the sum of powers phase region system 
can consist of different regions that have different slopes and/or different 
number of terms. For example, the linear phase region system, as shown in 
Figure 5A, instead of having a zero phase central region and two linearly 

1 0 changing regions, could have a zero phase central region and more than 

two linearly changing regions. Additional linear regions could direct the rays 
to different positions in order to better control misfocus, misfocus 
aberrations, and antialiasing characteristics. The size and shape of the 
misfocus PSFs can also be visualized and controlled by controlling regions of 

1 5 phase function and the corresponding rays. When only two asymmetric 

regions of the phase function are being used (as in Figure 5A), manual 
optimization is possible. When the number of regions exceeds two, then the 
number, size, and phase of each region is generally best determined by 
computer optimization. 

20 Figure 6 gives another view of the linear phase region and cubic phase 

region Wavefront Coding systems. The top graph of Figure 6 shows a 1D 
slice of phase functions describing the optics of both systems along one of 
the orthogonal axes. The 1 D form of this linear phase region system is: 

3.86 U(lxl/A x ) sign(x) [ (Ixl - A x )/(1- A x ) ] 

25 where 

A x = 2/3 
Ixl £ 1 

U(z) = 1 if z £ 1 , U(z) = 0 otherwise 
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sign(z) = +1 for z > 0, sign(z) = -1 otherwise 

The 1 D form for the cubic phase region system is specifically: 

5 U(lxl/A x ) sign(x) [ (Ixl -A x )/(1-A x ) ] 3 
A x = 1/3 

5 Ixl < 1 

U(z) - 1 if z < 1, U(z) = 0 otherwise 
sign(z) = +1 for z £ 0, sign(z) = -1 otherwise 

The 2D contours of constant phase for these phase functions are shown in 
the bottom of Figure 6. These contours clearly show that the phase within 
1 0 the central regions of these systems are constant; i.e. the central rays of 
the corresponding Wavefront Coding systems are not modified. The phase 
near the edges of the linear phase region system increases/decreases 
linearly while the phase near the edges of the cubic phase region system 
increases/decreases as a cubic function. 

1 5 Figure 7 shows the misfocus MTFs as a function of normalized misfocus for 

the system with no Wavefront Coding, the linear phase region system, and 
the cubic phase region system. The normalized misfocus values are y = { 0 
, 2 , 4 }, as used in Figure 3. The MTFs from both the linear phase region 
and cubic phase region systems are very close to each other (and thus are 

20 not specifically distinguished in Figure 7) and are very insensitive to misfocus, 
especially when compared to the system using no Wavefront Coding. 
Compare the heights of the improved Wavefront Coded MTFs of Figure 7 to 
those of the prior art cubic phase Wavefront Coded MTFs shown in Figure 3. 
The MTFs from both the linear phase region and the cubic phase region 

25 systems have higher MTFs than the conventional cubic phase function MTFs 
while all the Wavefront Coding MTFs are essentially insensitive to misfocus. 
In terms of antialiasing, the MTFs from the improved linear phase region and 
cubic phase region systems have higher values of MTFs, compared to the 
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prior art cubic phase MTFs, for the lower spatial frequencies that typically 
are not aliased, while also reducing the high spatial frequency optical power 
when compared to the system with no Wavefront Coding. 

Improved non-separable Wavefront Coded Optics can also be described. 
General non-separable Wavefront Coded Optics can be mathematically 
defined through the form: 

phase(p,8) = I Q(p/£^) G,( p,0 ) 

where 

Ipl < 1, 0 < 6 < 2 pi 
i = 1,2, ... N 
where Q(z) = 1 if z > 1, Q(z) = 0 otherwise 

0 < £^ < 1 

and where the sum is over the index i. The function Q(p/Q) is a zero/one 
function that allows the central rays from a region with radius less than Q 
to be unchanged. The function Gj( p,e) is a general phase function that has 
regions of positive and negative phase so that regions of rays are made to 
either cross the optical axis before the best focus image point, or after the 
best focus image point. 

One general example of improved non-separable Wavefront Coded Optics is 
mathematically described as : 

non-separable-phase(p,e) = X Q{p/QJ ol { p BI cos( w t 6 - ) 

where 

Ipl £ 1, 0 < 8 < 2 pi 

1 = 1,2,... N 

and where Q(z) = 1 If z * 1, Q(z) = 0 otherwise 
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o <o< < 1 

Another example of improved non-separable Wavefront Coded Optics is the 
non-separable sum of powers form given by: 

non-separable-sum-of-powers(p,8) = 
5 I Q(p/^)oC| sign angle (M0 - offset) [(p - 

where 

Ipl < 1, 0 < 6 < 2 pi 

1-1,2 N 

and where Q(z) = 1 if z > 1, Q(z) = 0 otherwise 
10 o <Q H <^ 

si 9 n angle( * ) = +1 if 0 S ° S P 1 ' Si 9 n angle( * ) = ' 1 Otherwise 

where the integer M controls the number of +/- sectors used and where the 
sum is over the variable i. 

An example of the improved non-separable sum of powers Wavefront Coded 

1 5 Optics is given in Figure 8. The upper left graph in Figure 8 is a contour plot 

of constant phase of the non-separable cubic phase sector system. The 
cubic phase sector system is the non-separable sum-of-powers system with 
one exponential term, with B - 3. M = 1 , and offset = pi/2. The 
parameters of this non-separable cubic phase sector system are: 

2 0 4.48 Q(p/&) sign angle (6 - pi/2) [(p - a)/(1-&)] 3 

fl SB 0.2 

Ipl < 1 , 0 < 6 < 2 pi 

and where 

Q(z) = 1 if z > 1 , Q(z) = 0 otherwise 
2 5 sign angle ( * ) = +1 if 0 < 0> < pi, sign angle ( * ) = -1 otherwise 

The top right graph of Figure 8 shows a 1D slice through the cubic phase 



1 7 



WO 03/021333 PCT/US02/09702 

sector system. The bottom graph of Figure 8 shows the misfocus MTFs of 
a system with no Wavefront Coding and with the non-separable cubic phase 
sector system. Again the normalized misfocus values are v = { 0 , 2 , 4 }, 
as used in Figures 3 and 7. These MTFs are 1D slices from MTFs formed 
with circular apertures. The MTFs for the system with the non-separable 
cubic phase sector system are seen to be very insensitive to misfocus 
effects, especially when compared to the MTFs from the system with no 
Wavefront Coding. 

Figures 9, 10, and 11 describe one example of using the improved 
Wavefront Coding Optics to control misfocus aberrations. This example 
shows control of field-dependent aberrations with a common Cooke triplet 
lens used for visible light. See Modem Optical Engineering (Warren J. Smith, 
McGraw-Hill, Inc, NY, 1990) for more information on the Cooke triplet. A 
drawing of this lens is given in Figure 9. The Noise Gain resulting from the 
improved Wavefront Coding Optics, as compared to prior art Wavefront 
Coding Optics, are shown to be a factor of 2.5 less than prior art systems 
with this lens. This reduced noise gain directly translates into final images 
that have 2.5 times less noise that prior art systems. 

This triplet uses all spherical surfaces except for the second surface of the 
second or middle lens, which contains the Wavefront Coding Optics. The 
prescription of the triplet lens without Wavefront Coding Optics is given by: 



Surface # Radius Thickness Glass Diameter 



Object 




Infinity 






1 


22.01359 


3.258956 


SK16 


17.10986 


2 


-435.7604 


6.007551 




15.98699 


3 


-22.21328 


0.9999746 




7.990108 


4* 


20.29192 


4.750409 




7.731661 


5 


79.6836 


2.952076 


SK16 


13.78242 
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6 


-18.3783 


42.18171 




14.5162 


Image 











Surface #4 is the aperture stop 
All dimensions are given in units of mm. The focal length of this lens is 
50mm, the full aperture F/# is 5, the half field of view is 20 degrees. This 
lens is used in a system with a digital gray scale detector. The pixel pitch is 
7.6 microns with 100% fill factor square pixels. The spatial frequency limit 
of this detector is 65.8 Ip/mm. The ambient temperature is considered 
fixed at 20 degrees C with one atmosphere of pressure. Ideal fabrication 
and assembly are also assumed. 

Even though the lens of Figure 9 is a fairly complicated multi-glass optical 
system, this system suffers from the monochromatic misfocus aberrations 
of spherical aberration, petzval curvature, and astigmatism that limit off- 
axis performance. If the lens were designed with a single optical material, 
instead of the two different materials actually used, the lens would also 
suffer increased chromatic aberration. If plastic optics were used in place 
of glass optics the lens would also suffer from temperature related 
misfocus effects. 

Figure 10A (prior art) shows the performance of the system of Figure 9 
operating at F/5 (wide open) and Figure 10B shows the system of Figure 9 
20 with the aperture stopped down to F/19. Each of the following graphs are 
related to measures of the lens with green light. A contour plot of the on- 
axis exit pupil optical path difference (OPD) for the F/5 system is given on 
the left of Figure 10A. The peak-to-valley OPD for this exit pupil has a value 
of approximately 0.7 X. The right plot of Figure 10A describes the MTFs for 
25 the F/5 system as a function of field angle for on-axis, 14 degrees, and 20 

degrees. Included in these MTFs, and all following MTFs, is the pixel MTF 

19 
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related to the 100% fill factor 7.6 micron square pixel. Only 1D slices along 
the horizontal axis of the 2D diffraction limited MTFs and the 2D Wavefront 
Coded MTFs are shown in Figures 10 and 11. The MTFs as a function of 
field angle are seen to vary greatly due to the presence of aberrations in 
5 the conventional F/5 system. 

The effect of stopping down the lens from F/5 to F/19 is seen in the graphs 
of Figure 10B. The peak-to-valley OPD of the on-axis stopped down exit 
pupil is reduced to 0.12 X. The MTFs as a function of field angle of the 
stopped down system are seen to change very little, compared to the full 
1 0 aperture system. Stopping down the aperture is one way of controlling the 

aberrations of this system, although at a reduction factor of (5/1 9) 2 or a 
loss of 93% of the optical power that was captured by the full aperture 
system. 

Figures 11A and 11B show the performance of the Cooke triplet of Figure 9, 
1 5 modified to utilize Wavefront Coding. The Wavefront Coding Optics for this 

example were modelled as being added to the aperture stop of the system 
which is at second surface of the second element, or surface #4. The 
surface height equations for the prior art Wavefront Coding cubic phase 
system used in Figure 11 A are: 

20 Z prior art (x,y) = 8.6061 E-5 {Ixl? + lyl 3 } 

Ixl < 5.45, lyl < 5.45 

where the surface height Z prjor art (x,y) is given in mm. The optical area that 

is used is a circle of radius 3.86 mm. A square aperture can also be used in 
practice. 

25 The surface equations for the improved Wavefront Coding cubic phase 

region system used in Figure 1 1 B are: 
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Zimproved< x .y) = 7.6535E-5 U(lxl) sign(x) [lxi-1] 3 
+ 7.653 E-5 U(lyl) sign(y) [lyl-1p 
Ixl < 5.45, lyl < 5.45 
U(z) = 1 if z > 1, U(z) = 0 otherwise 
5 sign(z) = +1 for z > 0, sign(z) = -1 otherwise 

where again the surface height Z lmproved (x,y) is given in mm, the optical area 
that is used is a circle of radius 3.86 mm, and a square aperture can also 
be used in practice. 

The graphs of Figure 11A describe use of the prior art rectangularly 
1 0 separable Wavefront Coding cubic phase Optics. The graphs of Figure 11 B 
describe use of the improved rectangularly separable cubic phase region 
Wavefront Coding Optics. Use of the improved Wavefront Coding Optics 
increases the height of the MTFs before image processing 110 and therefore 
drastically reduces the Noise Gain of the digital filters needed to give the 

1 5 ideal performance after processing, compared to the prior art optics. The 

improved MTFs also show increased antialiasing performance by having 
higher optical power below the detector spatial frequency cutoff, with 
greatly reduced MTFs beyond the detector cutoff, when compared to the 
conventional full aperture system. 

20 The prior art cubic phase system was designed to have a peak-to-valley 

OPD over the on-axis wide open (F/5) exit pupil of approximately 9.5 X. The 
resulting MTFs, before image processing, as a function of field angle are 
essentially constant, as shown in Figure 11 A, especially compared to the full 
aperture F/5 system without Wavefront Coding of Figure 10A. The 

2 5 diffraction limited MTF for this system was chosen as the desired system 

performance of the Wavefront Coding systems after Image Processing 110. 
In general, after image processing the resulting PSFs and MTFs in the 
Wavefront Coding system can have nearly any form. The Image Processing 
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function 110, for this example, implements a 2D linear filter to transform the 
Wavefront Coding MTFs before image processing to MTFs after image 
processing that closely match the ideal diffraction limited MTF below the 
detector spatial frequency cutoff of 65 Ip/mm. The resulting Noise Gain 
values of the digital filters are then used as a figure of merit to judge the 
two dimensional height of the MTFs compared to the desired diffraction- 
limited MTFs. For the prior art cubic phase system of Figure 11 A the Noise 
Gain of the resulting 2D digital filter is 8.1. 

The graphs of Figure 11 B describe use of the improved cubic phase region 
Wavefront Coding Optics. Since the stopped down F/19 system of Figure 
10B has suitable performance for rays within a 2 mm diameter aperture, 
the cubic phase region system is constant, or has zero phase, over the 
square aperture region of +/- 1.0 mm. This zero phase region corresponds 
to the aperture of the stopped down F/19 system. 

The zero phase region could just as easily be formed into a circle or other 
geometric shape depending on the application and processing to be used. 
The square aperture is more consistent with the rectangularly separable 
nature of the cubic phase region system than a non-separable circular 
region. 

The parameters of the cubic phase region system of Figure 1 1 B were 
designed to have similar insensitivity to off-axis misfocus aberrations as the 
prior art cubic phase system of Figure 11 A. This results in a peak-to-valley 
on-axis OPD also approximately 9.5 X. The contour plot of the exit pupil on 
the left of Figure 11B clearly shows a large region near the center of the 
aperture that has zero phase as compared to the contour plot from the 
prior art cubic phase system that is not optically constant over any region 
in the aperture. The MTFs, before image processing, as a function of field 
angle for the improved cubic phase region system are seen to be essentially 
constant. The height of the MTFs resulting from the improved cubic phase 
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region system is also much higher than those from the prior art cubic phase 
system below the spatial frequency band limit of the 7.6 micron detector, or 
65 Ip/mm. 

The Noise Gain of the 2D digital filter needed to match the performance of 
5 the MTFs after filtering to that of the diffraction-limited system has a value 

of approximately 3.2 for the improved cubic phase region system. Thus, the 
improved cubic phase region system produces nearly ideal performance in 
the control of the field dependent aberrations and also drastically reduces 
the digital filter Noise Gain from the prior art system by a factor of 
1 0 (8,1/3.2), or approximately 2.5. Thus, the additive noise power in the final 

images after Image Processing 110 will be 2.5 times larger with the prior 
art system than with the improved cubic phase region Wavefront Coding 
system. 

This large decrease in Noise Gain will result in substantially less noisy final 
1 5 images from the Improved Wavefront Coding Optics as compared to prior 

art optics. Or, for a constant amount of Noise Gain, the improved 
Wavefront Coding Optics can control a much larger degree of misfocus than 
can the prior art optics. 

Although not shown, this lens system when modified with Wavefront Coding 
20 also substantially corrects misfocus arising from chromatic effects and 

temperature related effects, as well as reducing the tolerance of the system 
to manufacture and assembly errors. 

What is claimed is: 
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CLAIMS 

1 . An improved Wavefront Coding system for imaging an object 
comprising: 

Wavefront Coding Optics (106, 601, 602) having an aperture and including - 

5 a central region, wherein the central region applies an essentially 

constant phase profile to light from the object passing through the 
central region, and 

a peripheral region disposed about the central region, wherein the 
peripheral region applies a phase profile to light from the object that 
1 o has alternating increasing and decreasing phase relative to the 

central region, 

wherein the phase profile applied by the Wavefront Coding Optics is 
aspheric and rotationally asymmetric, and 

wherein the phase profile applied by the Wavefront Coding Optics 
1 5 alters the optical transfer function of the Wavefront Coding system in 

such a way that the altered optical transfer function is substantially 
less sensitive to focus related aberrations than was the unaltered 
optical transfer function; 

a detector (108) for capturing an image from the Wavefront Coding Optics; 
20 and 

a post processing element (110) for processing the image captured by the 
detector by reversing the alteration to the optical transfer function of the 
Wavefront Coding system accomplished by the Wavefront Coding Optics. 

2. The Wavefront Coding system of claim 1, wherein the central region 
25 comprises a rectangular aperture, and the peripheral region forms a 

rectangular frame around the central region. 

3. The Wavefront Coding system of claim 1, wherein the central region 
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comprises a circular aperture, and the peripheral region forms a ring around 
the central region. 

4. The Wavefront Coding system of claim 1 , wherein the central region 
comprises a rectangular aperture, and the peripheral region forms a ring 

5 around the central region. 

5. The Wavefront Coding system of claim 1, wherein the peripheral 
region comprises more than one concentric zone. 

6. The Wavefront Coding system of claim 1, wherein the peripheral 
region of the Wavefront Coding Optics applies substantially linear 

1 0 functions. 

7. The Wavefront Coding system of claim 1, wherein the phase profile 
applied by the Wavefront Coding Optics substantially follows one of 
the following functions: 

(a) phase(x,y) = I [ U(lxl/A xl ) G xi (x) + U(lyl/A yi ) G yl (y) ] 

1 5 where 

1x1*1, lyl<1 
1 = 1,2, ... , N 
and where U(z) = 1 if z > 1, U(z) = 0 otherwise 
0<A xi <1, 0< A yi <1 

2 0 and where the sum is over the index i. 



(b) phase(x.y) = U(lxl/A x ) * sign(x)(lxl-A x )/(1-A x ) 

+ U(lyl/Ay) * sign(y)(lyl-A y )/(1-A y ) 

2 5 where 

Ixl < 1, lyl £ 1 
U(z) = 1 if z > 1, U(z) = 0 otherwise 
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sign{z) = +1 for z > 0, sign(z) = -1 otherwise 
0<A X <1 , 0< A y <1 . 



c) sum of powers phase(x.y) = 

I U(lxl/A xi ) Oj sign(x) [ (Ixl -A xi )/(1-A xl ) ] 0i + 
E U(lyl/A yi ) x , sign(y) [ (lyl -A yi )/(1-A yi ) )V 

where 

Ixl s 1, lyl £ 1. i = 1, 2 N 

U(z) = 1 if z > 1, U(z) = 0 otherwise, 
sign(z) = +1 for z 2 0, sign(z) = -1 otherwise, 
0<A xi <1, 0< Ayj^, and where the summations are over the index i. 



(d) phase(x,y) = V(\x\/\) a, sign(x) [ (Ixl -A x )/(1-A x ) ] 3 + 
U(lyl/A y ) x, sign(y) [ (lyl -A y )/(1-A y ) ] 3 

where 

Ixl s 1, lyl £ 1, 
U(z) = 1 if z £ 1, U(z) = 0 otherwise, 
sign(z) = +1 for z > 0, sign(z) = -1 otherwise, 
0<A X| <1, 0< A yi <1. 

(e) phase(p.e) = E Q(p/^) Gj( p,e ) 

where 

Ipl £ 1, 0 < 6 £ 2 pi 
i = 1, 2, ... N 
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where 

Q(z) = 1 if z > 1, Q(z) = 0 otherwise 
0 < < 1 

and where the sum is over the index i. 

5 

<f) phase(p.e) = X Q(p/Qj) a t p Bi cos( w s e - <!>, ). 

<g) phase(p.e) = X Q(p/£2i) sign angle ( M, 9 - offset ) G Pi (p) 

where sign angle ( * ) = +1 if 0 < <D< pi, sign aagle ( <D ) = -1 otherwise; 

1 0 where the integer M controls the number of +/- sectors used around the 

circle, and offset controls the rotation of the +/- sectors. 

(h) phase(p.e) = Z Qip/Clfa sign angle (M6 - offset) [(p - £2l)/(1 -ni)] 8i . 
1 5 (i) phase(p.e) = Q(p/Q) sign angle (M 6 - offset ) [(p - Q)/(1-£2)] 3 . 
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8. The method of designing improved Wavefront Coding systems 
comprising the steps of: 

designing Wavefront Coding Optics by - 

selecting a central region of the Wavefront Coding Optics such that 
the central region applies an essentially constant phase profile to light 
passing through it; 

selecting a peripheral region of the Wavefront Coding Optics such 
that the peripheral region applies a phase profile to light passing 
through it that alternately increases in phase and decreases in phase 
relative to the central region; 

arranging the peripheral region such that the Wavefront Coding 
Optics are aspherical and rotationally nonsymmetric; 

wherein the phase profile applied by the Wavefront Coding Optics 
alters the optical transfer function of the Wavefront Coding system in 
such a way that the altered optical transfer function is substantially 
less sensitive to focus related aberrations than was the unaltered 
optical transfer function; and 

selecting a post processing function for processing images from the 
Wavefront Coding Optics by reversing an alteration of the optical transfer 
function accomplished by the Wavefront Coding Optics. 

9. The method of claim 8 wherein the peripheral region applies 
substantially linear phase profiles. 

10. The method of claim 8 wherein the peripheral region applies 
substantially cubic phase profiles. 
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11. The method for increasing depth of field and controlling focus related 
aberrations in an imaging system for imaging an object at a detector, the 
method comprising the steps of: 

between the object and the detector, modifying the wavefront of light from 
5 the object; 

the wavefront modifying step including the steps of applying an essentially 
constant phase profile to light passing through a central region and applying 
a peripheral profile that alternately increases in phase and decreases in 
phase relative to the central region to light passing through a peripheral 
1 0 region disposed about the central region; 

the applying steps operating to result in an overall applied phase profile that 
is aspherical and rotationally nonsymmetric; 

the applying steps operating to result in an overall applied phase profile that 
alters the optical transfer function of the imaging system in such a way that 
1 5 the altered optical transfer function is substantially less sensitive to 
misfocus related aberrations than was the unaltered optical transfer 
function;and 

post processing the image captured by the detector by reversing the 
alteration of the optical transfer function accomplished by the the 
20 wavefront modifying step. 

12. The method of claim 11, wherein the step of applying a peripheral 
profile applies substantially linear functions. 

13. The method of claim 11, wherein the step of applying a peripheral 
profile applies substantially cubic functions. 
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Figure 2 (Prior Art) 
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Figure 5A 
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